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, Abstract 
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. We consider the grand canonical pressure for Coulombic matter with nuclear 

' charges ^ Z in a magnetic field B and at nonzero temperature. We prove that 

. its asymptotic limit as Z — > oo with B/Z^ — > can be obtained by minimizing 

' a Thomas-Fermi type pressure functional. 

^ : 1 Introduction 

This paper intends to add one more chapter to the saga of rigorous Thomas-Fermi 
theory in which Ehiott H. Lieb played a prominent role The issue is the derivation 
of Thomas Fermi theory at nonzero temperature in a strong magnetic field from 
quantum statistical mechanics. The asymptotic exactness of Thomas- Fermi theory 
I for Coulombic matter in its ground state was first proved by Lieb and Simon in the 

fundamental paper |2- A shorter proof, using coherent states, was given by Lieb in 
[2] , and several ideas in the present paper were inspired by that proof. 

Thomas- Fermi Theory for matter in extremely strong magnetic fields is important 
for the physics of neutron stars, cf. [H El El and references quoted therein. This 
theory was analyzed from the point of view of mathematical physics in (71 (HI El 
and its status as a limit of quantum mechanics in a certain parameter range firmly 
established; an extension of the asymptotics to inhomogeneous magnetic fields is in 
|inj . All these works are concerned with the ground state, but non-magnetic TF 
theory at temperatures T > has been treated in jlll 1121 [T^ . Magnetic TF theory 
at nonzero temperature was studied in J31 and applied to the equation of state for 
matter in the crust of a neutron star, but a proof of its asymptotic exactness has, 
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to the best of our knowledge, not been published until now. The proof we give here 
brings together techniques from jlll I12j and [Tj |Sj with several amendments and 
additions. 

We start our discussion with some heuristic considerations. A possible point 
of departure for a motivation of TF theory, both at T = and T > 0, is the 
thermodynamic relation between the particle density p, the chemical potential //, 
and the pressure P for a homogeneous gas of noninteracting particles (electrons): 

/, = ap(M)/a/z=:P'(/z). (1.1) 

The next step is to consider electrons that interact with with each other via Coulomb 
forces and also with an external potential V (arising from nuclei in fixed positions as 
well as a confining potential that prevents the electrons from escaping to infinity). 
The electron density now depends on the position x. The TF theory is formally 
obtained from (ll.lj) by replacing the constant density p by a position dependent den- 
sity /o(x) > and the chemical potential by a position dependent chemical potential 
//(x), imposing as an equilibrium condition that the total electrochemical potential 

= /i(x) + Fp(x), (1.2) 

with 

yp(x) = y(x) + p*|xr\ (1.3) 

should be independent of x. The result is the Thomas-Fermi equation: 

/.(x) = P'(pTF-V^p(x)). (1.4) 

For given P, V and /xtf this is a nonlinear integral equation for yo(x). The total 
particle number is 

j p{x)dx = N. (1.5) 

The equation (|1.4jl is the variational equation associated with the minimization 
problem for the TF pressure functional of the density,* 

V[p] = j P {fiTF - I^p(x)) rfx + D{p, p), (1.6) 

with 

The minimum of H1.6() over all nonnegative functions p will be called the TF pressure. 

*Instead of considering P as a functional of the density, it could equivalently be considered 
as a functional of the potential Vp. Note that p and hence D{p, p) is determined by Vp because 
47rp(x) — — Ap * |x|~^ — A(V — Vp(x)). While this point of view (which is related to that of Firsov 
|15| in standard TF theory) may be more natural if V is regarded as a Legendre transformation of 
the free energy functional J- ILIOL we find it convenient in the present context to regard "P as a 
functional of p. 
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An alternative form of the TF equation is obtained if one replaces P by its 
Legendre transform, the free energy density 

f{p) = snp{fip-P{fi)}. (1.8) 

From (|1.8|) it follows that df /dp =: f is the inverse of P'. Hence the TF equation 
H1.4|) can also be written 

/'(p(x))+yp(x))=/ZTF. (1.9) 

This form of the TF equation is associated with the minimization problem for the 
free energy functional 

Hp] = I {/(P(x)) + ■t^(x))p(x)} dx + D{p, p) (1.10) 

with (|1.5() as a subsidiary condition and ;Utf as a Lagrange multiplier. 

In the sequel we shall investigate the functional (|1.6|) . Our main result is its 
asymptotic exactness in the case of electrons in a constant magnetic field and at 
nonzero temperature. This amounts to taking a semiclassical limit of the grand 
canonical pressure for the quantum mechanical many-body problem. The corre- 
sponding problem for the free energy, i.e., the canonical partition function, is tech- 
nically more involved and will not be tackled here. 

We now introduce some notation that will allow us to state the results precisely. 

The many-body Hamiltonian considered in this paper is 

N 

i^w,z,B = X]{[pi+A(^i))-'^i]'+^^.s(^*)}+ ixi-xj-r' (1.11) 

i=l l<i<j<N 

Here p = — iV, A(x) = i(— i?X2, 0) is the vector potential of a homogeneous 
magnetic field of strength B in the xa-direction, and a- is the vector of Pauli matrices. 
Atomic units are chosen so that h = 2m = e = 1, and the temperature unit is such 
that the Boltzmann constant k is also 1. The external potential is 

K 

Vz,b{?^) = -^E |x-k I + ^^"'^(^"'^) (1-12) 
k=i ''^ 

where is a confining potential that tends sufficiently rapidly to oo for |x| ^ cx) 
so that exp(— I^(x)/T) is integrable for all T > 0. It will also be assumed to satisfy 
some regularity conditions stated later. The are fixed positions of nuclei with 
fixed charges Zk <: 1 which are scaled by an overall parameter Z. The length scaling 
factor 

£ = £z,B = Z-^l\\ + (5/Z4/3)]-2/5 (1.13) 
is the one appropriate for TF atoms in a magnetic field, cf. |S]. 



4 



The Hamiltonian Hn,z,b operates on the A^-electron Hilbert space of antisym- 
metric wave functions in space and spin variables: 

Hn = A^L\R^,C^). (1.14) 

The corresponding Fock space is 

oo 

n = ^nN (1.15) 

N=0 

with TCo = C. If An are operators on TCn, N = 0,1, . . . , we denote the operator 



^qAn on H by A. In particular, 



Hz,B = Hn,z,b. (1.16) 

7V=0 

The grand canonical pressure at chemical potential fi and temperature T is* 

P^^{fi,B,T,Z) = T Intr exp[-{Hz,B - f^N)/T]. (1.17) 
The free 1-particle Hamiltonian 

Ho = [(p + A(x)) . (1.18) 

has the Landau spectrum 

e^ip) = 2Bu + p"^ , u = 0,l,..., peR (1.19) 
with degeneracy (pro unit area in the (xi, X2)-plane) 



duiB) 



B/{2tt) if iy = 
B/iT if u>l. 



The pressure of a free electron gas at temperature T in a magnetic field and with 
chemical potential fi is 

Pt,b(^) = tY,MB) / ln[l + exp{-(e,(p) - fi)/T}] dp. (1.20) 

The magnetic Thomas-Fermi pressure functional is obtained by taking P in 1)1.6(1 to 
be (tr^ : 

7'^T''[p;Z,/i,r,S] = j Pt,b{i^-Vz,b,M)) d^^ + D{p,p) (1.21) 

with 

Vz,bA^) = VzM^) + P * l^r'. (1.22) 
The pressure according to MTF theory is 

P^'^^{Z,n,T,B) := miV^^^[p;fi,T,B,Z]. (1.23) 

p>0 

We can now state our main result: 



*Note that n denotes here the total electrochemical potential and not the chemical potential of 
the free electron gas denoted by the same letter in in jLlJ and lfL20jl . 
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Theorem 1.1 (The MTF pressure is a hmit of the QM pressure). // 

Z,fi,T oo with 11/ {Zi'^^) and T / {Zl'j^^) fixed hut B/Z^ — > 0, where Ib,z is 
given by U.lcl\) . then 

P^^(n,T,B,Z) , , 



P^^^{li,T,B,Z) 



1 



The main steps in the proof are as follows. In the next section we discuss the MTF 
functional in more detail, in particular the existence and uniqueness of a minimizer 
and its properties. In Section 3 we consider a Hamiltonian with a mean field and the 
corresponding pressure functional. Section 4 contains the basic semiclassical limit 
theorem, which is proved using the magnetic coherent states introduced in 7, and 
[H] . As a corollary one obtains the convergence of the mean field pressure to the MTF 
pressure. Here it is important that the error bounds in the semiclassical theorem are 
uniform in the relevant parameters. In the last section the proof of Theorem 1.1 is 
completed by estimating the many-body pressure in terms of the mean field pressure. 



2 The MTF pressure functional 

In this section we collect some basic facts about the pressure functional ()1.21|) . We 



start with some formulas and estimates for the pressure (|1.20j) of the free electron 
gas in a magnetic field. It can also be written as 



G{e) 



/o exp{(e-/i)/r} + l 
with the integrated density of states 

B 



de 



G{e) 



21/2^2 



(.1/2 + 2^ |e - 2Bv 



1/2 



v=0 



(2.1) 



(2.2) 



Here |s|+ := \s\ for s > and otherwise. The corresponding formulas for Pqp ^ 
dPT,B/d^ are 



00 „ 



1 



i/=0 

00 



exp{(e,(p)-/i)/T} + l 



dp 



/o exp{(e - /x)/r} + 1 
The scaling properties of Pt,b can be seen by writing 



de. 



(2.3) 
(2.4) 



and 



PtAi^) 



2i/2vr2 

Bpy^ 



h/2 it^/T) + 2 A/2 (/i - {2Bi^/T)) 



1-1/2 it^/T) + 2 J2 1-1/2 (/« - {■^Bv/T)) 



(2.5) 



(2.6) 
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with 

/■oo k 

Using the inequahties s) < (e^ + < 0{—s) + e~l'*l) we obtain from (|2.1() the 
following simple estimates, treating the sum over > 1 as a Riemannian approxima- 
tion for an integral: 

c(i3|^|f + |^|f ) <Pt,b(m) 

< C (slMlf + \p\T + e-\>^\/^{BT^/^ + r5/2)) (2.8) 

with constants c > and C < oo. In the same way we obtain 
c'[B\f,\'f + \f,\f)<P/r,M 

< C {B\^l\^_l^ + 1^15./' + e-i'^i/^(i?ri/2 + r3/2)) . (2.9) 

As domain of definition for the magnetic pressure functional we take 

M = {p: p(x) > 0, D{p, p) < cx)}. (2.10) 

To see that V^'^^[p; p,T, B , Z] is well defined and < oo for /) G we use (ITKll 
together with 

Lemma 2.1 (Coulomb bound). If p > with D{p,p) < oo, then the potential 
fp(x) = p* [xj^-^ is in L^(M^); in fact 

\\vp\\l<{const.)D{p,p). (2.11) 

Proof. Since 47r/9 = — V^Up, we can write D{p,p) = (const.) |[Vt>p||2- We then use the 
Sobolev inequality in [TB], i.e., ||/||6 < (const.) ||V/||2. □ 

Prom this lemma follows that p — Vz,B,p £ -^^''^(I^^)ioc ^-nd hence also in 
L^/'^{'R%oc for p e M. Moreover, since Ty(x) tends to oo if x ^ oo while the 
negative potential from the nuclei tends to zero, there is an ii > such that 
Ia* - yz,B,p\+ = and \p - Vz,B,p\ > Z£-'^W{e-'^x) - (const.) for |x| > R. The 
finiteness of V^'^^lp; p,T, B, Z] now follows from H2.8() by splitting the integration 
domain in ()1.21() into |x| < i? and |x| > R. 

The MTF pressure functional is nonnegative and strictly convex since D{p,p) is 
strictly convex, p Vp is linear and s i— > ln(l + exp(— s)) strictly convex. Moreover, 
if Ai is equipped with the topology defined by the Hilbert norm D{p,p)'^/^, then 
•pMTF jg -^^gakly lower semicontinuous. This can be seen by representing the convex 
functional as a supremum of affine, weakly continuous functionals in a similar way as 
in ^5], (4.1.10). This continuity and strict convexity implies by standard arguments 
(c.f., e.g. p]) that the functional has a unique minimizer f^"^^ g Jv[ (depending on 
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the parameters fj,,T, B, Z). It is the unique solution to the MTF equation, i.e., the 
variational equation for the minimization problem, 

p(x) = - Vz,bA^)). (2.12) 

Eq. (|2.5)) implies the following scaling property of the MTF functional: 

V^^^[p;fi,T,B,Z] = z2rip^T^[/5;A,T,/3] (2.13) 

where 

P = B/Z^/^, ^ = Z-i/3(l+/3)-2/5, (2.14) 

/9(x) = zr3p(r^x), (2.15) 
fi = fi/izr'), f = T/{zr') (2.16) 

and 

P^^TF[~. /?] = (1 + /J)-3/5 y _ V>(x))dx + i?(p, p) (2.17) 

with 

B = S/(Zr^) =/?(l + /?)-2/5^ (2.18) 

and 

V>(x) = -$c(x) + Ty(x) +p* |xr^ (2.19) 

with 

<^>c(x)=X;^^- (2-20) 
X — Xi. 

fc=i ' 

We can also include the limiting cases /? = and (3 = oo: The case /3 = is 
just the temperature dependent TF without magnetic field considered in [1111121 [T^ . 
while P ^ oo means that only the lowest Landau level contributes. It is the T 
dependent version of the STF theory in 8 . In fact, as is easily seen from 1)2. 5() . 

f3/2 

hm (1 + /3)-3/5p^ ^(p) = ^^/,/2(^/r) =: P^ip) (2.21) 

and hence 

lim V'^^''[p;p,f,P]= [ P:^{p-Vp{^))d^ + D{p,p)=:V^[p-p,f] (2.22) 

As last topic in this section we derive some uniform bounds for the minimizing 
density of the MTF functional and the corresponding Coulomb potential. These 
bounds will be needed for the semiclassical limit theorem in Section 4. 
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Lemma 2.2 (Bounds for MTF minimizer). Let pp he the minimizer of p ^ 
V^'^^ [p; fl, f, P] for fixed jl, f. Then 

(i) The potentials W/3(x) = pp* |x|~^ are bounded in L^(M) uniformly in (3. 

(a) pp € L*'(M^) for 1 < p < 2 and \\pp\\p is uniformly bounded in 0. 

(Hi) Let jr = r~^j{-K/r) where < j € C^(M^) satisfies J j(x)(ix = 1. Then for all 

R < oo and p < 3 

/ |w/3(x) - * |x|-Ydx ^ (2.23) 

J|x|<_R 

uniformly in (3, as r ^ 0. 

Proof, (i) This follows from Lemma l2. II and the bound 

D{pprpp)<V^^^W,f^,f,(i]<V^^^[Q-f,,f,P]. (2.24) 

The right side is continuous in /? > and converges to 7'°°[0; /i, T] < oo for /? — > cx). 
Hence D{pp,pp) is uniformly bounded in (3. 

(ii) The MTF equation for the scaled density is 

p^(x) = (1 + /3)-3/5p^_^(/i + cl>^(x) - W{^) - pp * Ixr^). (2.25) 

Now, since P' is monotonously increasing, the right side is bounded by 

(1+ P)-^/^pL^^{p + ^c{^) - W{^)) 

< C (lA + $c(x) - W{-x))tl^ + exp(-|/i + $c(x) - H^(x))|/r)) (2.26) 

which is in LP{M?) for ah 1 < p < 2. □ 

(iii) This is proved in the same way as Proposition 4.19 (iii) in 8 , using Jensen's 
and Young's inequalities together with (i) and (ii). Note that the Coulomb potential 
|x|-i is in LP(M3)ioc for p < 3. 

3 Mean field theory 

For p E we define a mean field Hamiltonian by 

Hz,B,p = [(P + A(x)) • + Vz,bA^) (3-1) 
and a mean field pressure functional by 

P-f[p;/i,r,S,Z] =rtrln[l + exp{-(Fz,B,p-^)/r}]+D(p,p). (3.2) 
Note that the first term is equal to 

T lntrexp[-(Fz,ij,p - pN)/T] (3.3) 
where Hz,B,p is the second quantization of Hz,B,p- 
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By exactly the same methods as in |12j . (4.1.10)-(4.1.13), one can show that ()3.2() 
is strictly convex and weakly lower semicontinuous on A4 and that the minimizer, 
p™^, is the unique solution of the self-consistent (Hartree) equation 

p(x) = 2 (x {exp{{Hz,B,p - f^)/T} + ly' x\ . (3.4) 



The right side is here the diagonal of the integral kernel of the trace class opertor 
{exp{{Hz^B,p — y)/T} + Next we introduce the unitary operator 

([/,V) (x) = r3/V(rix) (3.5) 

for il) G L^(M^;C^) with £ given by (|1.13l) . It transforms the Hamiltonian according 
to 

Ui^Hz,B,pUe = {Zr^)Hh,b,p (3.6) 

with 

Hh,b,p = [{hp + foa(x)) • + V>(x). (3.7) 
Here a(x) = i(— X2,xi,0), Vp is defined by ()2.19|) . and 

h = r^/2^-^/2 = Z-l/='(l+/?)l/^ (3.8) 
b = Bf/^Z'^^^ = Z^/^P{l+P)-^/\ (3.9) 

Since the trace is invariant under a unitary transformation we obtain 

T tr ln(l + exp{-{Hz,B,p - l^)/T}) = (3.10) 
(Zr^) f tr ln(l + exY>{-{Hh,b,p " ~^)/T}), (3.11) 

where T and fl are given by (|2.16jl . Hence, 

V'^'Ip; II, T, B, Z] = Zr^ftv ln(l + eM-{Hh,b,P " f^)/T}) + Z^r^D{p, p) 

Z-^ftr ln(l + exp{-{Hh,b,p - p)/f}) + D{p, p)\{2,.l2) 



z'^r^ 



In the next section we study the semiclassical limit — > of (|3.12() . which is equiva- 
lent to Z ^ oo with B/Z'^ 0. The minimizer depends on b and h, and in order to 
obtain a semiclassical limit for the minima we need uniform bounds on the minimiz- 
ers and the corresponding Coulomb potentials. The following lemma is analogous to 
Lemma 12.21 

Lemma 3.1 (Bounds for mean field minimizer). Let pb,h be the minimizer, for 
fixed p, T, of 

V'^'[p;p,f,b,h] = [z-'fti ln{l + exp{-{Hh,b,p-p)/f}) + D{p,p)\ . (3.13) 
Then 

(i) The potentials Vb^h{^) = Pb,h * |x|^^ are hounded in L^(IR) uniformly in b,h. 
(a) Pb,h S -L^(M^) and ||p6,7i||i is uniformly bounded in b,h. 
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(in) Let jV = r ^ii^/f) where < j € C^(M^) satisfies J j(x)(ix = 1. Then for all 
R < oo and p < 3 

[ \vb,h{^) - Vb,h * [xj-i Pdx ^ (3.14) 

J|x|<i? 

uniformly in b, h, as r ^ 0. 

Proof, (i) As in Lemma 12.21 the essential point is that D{pb^h,Pb.h) is uniformly 
bounded. In fact, 

D{pb,h,Pb,h) < r'^'[pb,h;p,f,p] < P"^f[0;A,r,/3]. (3.15) 

From the magnetic Lieb-Thirring inequality ^8- it follows in the same way as in |12j . 
(4.1,47), that the pressure "P™^ is, up to a constant factor, bounded by the MTF 
pressure. The uniform bound thus follows in the same way as in Lemma 12.21 (i). 
(ii) The Hartree equation for pb.h is 



Pb,h{^) = 2Z-U^ 



x). (3.16) 



Since p * |x|~^ > 0, the integral of the right side is bounded by the trace of the 

operator 2Z^^ (^exp{{Hh,b,Q — A)/^'} + 1^ i which can, again by a magnetic Lieb- 
Thirring inequality, be bounded by the corresponding semiclassical expression, i.e., 
(1 + J {P'^ j^{p + ^^(x) — VF(x))(ix. This is bounded uniformly in the param- 

eters, by the same estimate as in Lemma 12.21 fii). 

(iii) This is proved in the same way as in Proposition 4.19 in |H], using (i) and (ii). □ 



4 Semiclassics 

We consider generally the operator 

Hh,b,v = [(/ip + ^'a(x)) • af + v{-k) (4.1) 
with /i > 0, 6 G M, a(x) = i(— X2, a^i, 0), and 

7;(x) = 't;i(x) V2(x) (4.2) 

where vi € L^f^^ and V2 is continuous with V2{^) — > oo for [x[ oo. We also impose 
some further conditions on vi and V2 that are described below. In the application to 
(PTTTl) . 

z;(x) = y^(x)-/i, (4.3) 

and we shall take 

t>i(x) = -^'c(x) -Fp * jxp-^ (4.4) 

and 

V2{^) = W{^)-p. (4.5) 
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The goal is to find an asymptotic approximation for 

P'^{h,b,v,T) ■.= TtTln{l + exp{-Hl,^Jr)) (4.6) 

as /i — * 0. 

The specific conditions for vi and V2 are: 

1. We assume that outside some compact baU, Bn^^ = {x : [x] < Rq}, vi is 
subharmonic and the negative part, |wi(x)|_ := |t'i(x)| for vi(x) < and 
otherwise, is continuous and tends to zero at oo. This is fulfilled for (|4.4|) if Rq > 
maxfc |Xk|. Subharmonicity implies the following property that is convenient 
for the proof of the upper bound: Let g G C^(M^) with f (jr^(x)(ix = 1 and 
define 5r(x) = r~^/'^g{x/r) for r > 0. Then 

vi{^) - VI * g^.{^) > (4.7) 

for X ^ and r small enough. 



2. We assume that V2 tends to oo sufficiently rapidly, so that e '"^^'^/'^ g 
for all T > 0. (This implies in particular that e~'"^^'^^'^ G LP(M^) for all p > 1.) 

3. We assume that V2 is sufficiently regular so that following holds: Let g G 
C^(R3) with / c/2(x)(ix = 1 and define ^^(x) = r-3/2^(x/r) for r > 0. Then 
we assume that there exists a continuous function V2 such that 

V2*9r= V2- (4.8) 

Moreover, limr_^o^2(^) ~ V2{'x) for all x, and e~^2( )/'^ < f(-,T) G L^(M^) for r 
sufficiently small. The same conditions should be fulfilled for V2,r '■= V2* g1- 

These conditions on V2 are not very restrictive. In fact, Eq. 1)4. 8() has a solution 
for all C°° functions V2 by Sect. 16.5, and by Fourier transform it easy to 
check explicitly that the conditions hold, e.g., for all polynomials. 

With Pt,_b(/u) the pressure of the free electron gas, cf. ()1.2U|) and 1)2. 5(1 . we define 

P{w-h,h,T) ■.= h'''Pr,hb{w) (4.9) 

f'OO 

= T^dy{h,b) / In [1 + exp{-((/ip)2 + 2i//i5- 'u;)/r}] dp(4.10) 

n J —00 



and 



P'"'\h,h,v,T) := j P{-vix);h,b,T)dx. (4.11) 

Then the following holds: 
Theorem 4.1 (Semiclassical limit theorem). For fixed t and v 



P^(h,b,v,T) 

lim ; ' ; = 1 (4.12) 



uniformly in b. 



12 



(4.13) 
(4.14) 



Proof. We shall make use of convexity of the function 

(/)(s) =rln(l + e-^/"), 

which implies the inequality 

(f){s + t) > ^{s) + (t>'{s)t. 
Moreover, if L, M and L + M are self-adjoint operators then for all a > 1 

tr 4>{L + M) < tr (/)(L) + - [tr 0(L + aM) - tr0(L)] . (4.15) 

a 

Another important tool for the proof are the magnetic coherent operators 
n(i^, u,p) with V = 0,1,2, .. ., u € and p G M introduced in [Jj and [H]. These 
operators fulfill the following conditions, cf. Eqs. (3.16)-(3.23) in "W: 

(4.16) 
(4.17) 

(4.18) 
(4.19) 



n{u,u,p) > 0. 

^ / / Il{i^,u,p)dpdu = 1, 

V 

trn(z^, u,p) = dy{h,h). 



iiHh^i,^Ji{v,\i,p) = dy{h,b) ep^y{h,b) +v * gl{vL) + —Ig 



with £p,y{h, h) = {hpY + 2bhv and Ig = f (V(7)^(x) dx. Moreover, 



h,b,v 



^p,i^ih, b) +vi{u) + V2iu) -^Ig n{u,u,p)dudp 
+ (^;i -vi* gl). 



(4.20) 

We also use the following inequalities that are easy consequences of (|4.17|) - ()4.18() 
and convexity of (p: 

tT^{Hh,b,v)>Y, j j du{h,b)(t>(^-^^-^tT{Hh,b,vU{u,u,p))^ dudp; (4.21) 

^/ J ^^^^^''^'P^^ ^^'^P-^'^^i^ J y /(^'"'P)n(i^,u,P)c^U(ip^ , (4 

where Eq. 1)4. 22() holds for all sufficiently regular functions / such that both sides 
well defined. 

The last preparatory step is to note that Eqs. 1)2. 8|) and ()2.9|) imply the following 
estimates for P{w;h,p,T): 

c{h'%\wtl'^ + h-'^\wtl'^^ < P{w;h,p,T) 

< C {h-H\w\%^ + h-^\wtl^ + e-\^\l^{h-%T^/^ + /i-3t5/2)) (4.23) 
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and 

P'{w; h,p,T) < C (h-H\wtl'^ + h-^wtl"^ + e-^'^^l^ih-HT^I^ + h-^T^^/^)) . (4.24) 

Eq. H4.23() implies in particular that for fixed r and v 

P'"'\h,h,v,T) ^ {h~H + h-^). (4.25) 

In order to prove (|4.12|) we thus have to estimate r In (1 + exp(— if/j^^ t,/r)) from 
above and below by P'*'^'(/i, 5, u, r) with errors that are small compared to (|4.25p . 

4.1 Lower bound 

We use (|4.19|) and begin by writing 



/j2 

£p^y{h,b) + v*gl{\i) + -j^Ig = A{v,\i,p) + B{i/,u,p) (4.26) 



with 



A{iy,u,p) = ep^^{h,b) +vi{u) +V2* gliu), 
B{u,u,p) = —Ig + {vi*g'^{u)-vi{u)). 



According to (|4.2H) and (|4.14|) we have 

P^{h,b,v,T) > [ d^ih,b)(p{A + B)dudp 

u 

-5^/ j du{Kb){(P{A) + (P'{A)B] dMdp. 

u 

We have to estimate the last term, i.e., the integral over (j)'{A)B. We use 1)4. 24() and 
the assumption that |t;i(x)|_ — > and V2{x) ^ oo if |x| oo, so | — u(x)|+ = 
for X outside some ball Br = {x : |x| < R}. Since V2 is continuous and V2 * g^ 
therefore bounded on Bpt, uniformly in r, the terms involving Ig are, after division 
by {h^'^b + h^^), bounded by 

(const.)!/ |7;i(x)p/2^x + l+ / e^p{-V2 * gr{^)/T)dy:] /r^. (4.27) 

\J\yc\<R J\x\>R J 

Since vi € L^^^ C L'^/^^ and exp(— ^2 * dr/''') is bounded by an function, inde- 
pendent of r, ()4.27p tends to with h if r = , < 6 < 1. The terms involving 
vi * gr — vi are, up to the factor {h^'^b + /i~'^), bounded by 



(const.) / {\vi{x)\^/'^dx + l)\\vi* g^{x) -vi{x)\dx< 

J\x\<R 

(const.) ( [ (|fi(x)[^/2^x + l)|dx) ( [ \vi* gf.{x) - vi{x)f^^d^ 



(4.28) 



5/2 

which tends to zero with r because vi G Lj^^ . 



14 



4.2 Upper bound 

Here we use (|4.20p , (|4.15p and (|4.22|) . In addition we need the Lieb-Thirring inequal- 
ity for a constant magnetic field (see [Sj), from which it follows in the same way as 
in [H], Ex. 1 (4.1, 47), that 

pQ(/i, b, V, t) < {const. )P''''\h, b, v, r). (4.29) 

Now by (|On|) we can write Hl^^ = L + M with 



Sp^„{h, b) + i;i(u) + V2{\i) - -^Ig 



n(i/, u,p) du dp 



and 

M = {vi- vi *gl). 
According to (|4.15|) we thus have 

P^{h,b,v,T) < ti (f){L) + - [tr (t>{L + aM) - tr(t>{L)] (4.30) 

a 

for all a > 1, and by (|4.22|) we have 

tr <P{L) < P''''\h, b, vi+vl^- h^Ig/r^, r). (4.31) 

If r = , < 6 < 1, then it follows from the properties of V2 and the dominated 
convergence theorem that the right side of 1)4. 31(1 converges to P^''\h,b,v,T) (in the 
sense that the ratio tends to 1), if /i ^ 0. 

We thus have to show that it is possible to let a — cxd as r — s- 0, in such a way 
that tr(/)(L + aM) stays bounded by (const.) • {h~^ + 6/i"^). 

Now L + aM = Hl f^^^ + {a - 1)L and thus 

tr (p{L + aM) = P^{h,b,v + (a - l)[vi - vi * g2],T). (4.32) 
To estimate this we use the inequality (|4.29|) which gives 

tr 0(L + aM) < (const.) j P'"'\h,b, -v{x) - {a-l)[vi{x) -vi* g'^{x)],T) dx. (4.33) 

We estimate this further using (|4.23|) as well as the assumptions on vi and V2, which 
imply that, outside some compact ball Bji, 

I - v{x) - (a - l)[t;i(x) - VI * g^{x)\+ = (4.34) 

and 

I — f (x) — (a — l)[t;i(x) — vi * gr{^)\ > \v2{x)\ — (const.). (4.35) 

Now 

/ [i;i(x) -vi* gl{x)f/^dx < ip{r) (4.36) 

J|x|<_R 

for some function ip with ip{r) — > if r — > 0. If we choose a = 1 + ip{r)~^ it follows 
from 1)4. 23() that tr ip(L + aM) is bounded by (const.) • {h~^ + bh~'^) for all sufficiently 
small r. □ 
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In the course of the proof we have shown that (for fixed r and v) 



P^{h,b,V,T) - P'^\h,b, V, T) 
and hence, for h and b defined in H3.8() . H3.9|l . 



h^P^ih,b,V,T) 



-Pr,/i6(-^^(x)(ix 



0(1 + hb) = oi{l + p) 



3/5 N 



(4.37) 



(4.38) 



Note also that = Z-^{1 + pf/^, hb = /?(1 + /3)-2/5 = b. We now apply this to 
compare the MFT functional (|1.2H) with the mean field functional (|3.2|) . making use 
of the scalings 1)2. 13(1 and ()3.12|) . At this point we need to assume that the confining 
potential W satisfies the regularity conditions stated for V2 above. Without further 
ado we obtain 



Corollary 4.1. [Convergence of functionals] If Z 
and T are fixed, then 



^ (X) while B/Z^ and p, fi 
Z-^lV'-'[p-f,,T,B,Z]- P^^TF ip.f,^f,(3] 



0. 



(4.39) 



For the application to Theorem 1.1, however, we need more than convergence of 
the functionals for fixed p, namely the convergence of the minima:* 

Corollary 4.2. [Convergence of minima] If Z —* oo while B/Z^ — > and p, T are 
fixed, then 

(4.40) 



Z-'H P'^^ifi, T,B,Z)- P^^^^^ (/i, T, (3) 



MTF/r, 



0. 



Proof. Let pfs be the minimizer of p i— > T^'^^lp; fl,f, P] for fixed p, f. Let pp be 
the corresponding unsealed density given by (|2.15|) . If Z ^ oo with /? fixed (which 
implies B/Z^ 0) it follows from Corollary 14. II that 



Z-^iP'^\fi,T,B,Z) < Z-^eV'^'[pp;p,T,B,Z] 

- V''^''[pf,;p,f,P] = P^^^{f^,f,P). (4.41) 

The general condition B jZ'^ 0, however, allows (3 and hence pf^ to vary as Z ^ oo 
{h 0) so one must check that the error terms in the semiclassical proof are uniform 
in (5. These terms involve J[xi<ij l^i(^) — ^^i * fi'r (x)!^''^'^^ with vi = — $c + Pf3*\'\~^ 
and the required uniformity follows from Lemma 12.21 The converse inequality. 



\imZ-'HP'"\p,T,B,Z) > P'^''''{fi,T,P) 



MTFt 



(4.42) 



follows in the same way by noting that the minimizer ph^b of the mean field functional 
enters in in the error terms of the upper bound only through L^-norms oi Vhh = 



Ph,b 



1-1 



and Vhi 



Vh,b * 



g'l that are uniformly bounded by Lemma l3. II 



□ 



*Pointwise convergence of functionals does in general not imply convergence of their infima, even 
if the functionals are strictly convex and their minimizers lie in a compact set. This can be seen 
from the following example : Take J-n{x) — x„ + ||a;|p on the Hilbert space P of square summable 
sequences x — {xi,X2, . . . ). Then J-„{x) — > —: J-{x) for all x and the infimum of T is zero. 

The infimum of !F„, on the other hand, is —1/4. All infima are attained in the weakly compact unit 
ball in f. 
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5 Proof of the QM limit theorem 

In this section we complete the proof of Theorem 1.1. We use the notation explained 
in the Introduction, in particular Hz.b for the Hamiltonian ()1.16|) on Fock space and 
Hz,B,p for the second quantization of the mean field Hamiltonian ()3.ip . 



5.1 Upper bound 

We use the inequality ^Sj 

N , N 



E^^^E/ T^d^-D{p,p)-3,m^N-^ [ p'/\^)d^ (5.1) 
TT^I^i-Xjl Z~^^M3 |x-Xi| 57 71.3 



t<j •' 1=1 



that holds for all 7 > and all p e M n (M.^) n L^/^{R^). Eq. ((SU) imphes 

Hz,B > Hz,B,p - 3.687iV - D{p, p) - C^^p (5.2) 
with C^^p = ^ Jjg3 /9^/^(x) dx. Since A > B implies tr < tr , we get 



upper bound on the grand canonical pressure: 

P^^{p,T,B,Z) < Tin (ti e'^^^'^'^-'"'^^/^) +D{p,p)+C^,p 



an 



= Ttr ln(l + e-(^z.s,p-/.^)/T) ^ ^^^^ ^ (j^^^ (5 3) 
= V''''[p;p^,T,B,Z] + C^,p (5.4) 

where p^ := p + 8.687. We now apply the unitary scaling (|3.6() explained in Section 
3 and note that 

ll^ = jl + (Z~^^)3.687 (5.5) 

and 

C,,p = {ZH-')^-I^Jp'/H^)d.. (5.6) 

We choose 7 as a function of Z and B such that Z~^£^ — > but Z^^^i^ — > 00 if 
Z —* 00, which is fulfilled as long as 

[1 + (B/Z4/3)]2/5 « J « Z*/^[l + (5/Z4/3)]2/5. (5.7) 

Then, according to (|3.12|) and Corollary 14.11 we get as Z — > 00 with B/Z^ and 
p, T fixed: 

llm{Z-^l)P^^{p, T, B, Z) < V^^^[p; p, f, /3] (5.8) 

for each /5 € A1 H L^(M^) n L^^^{M^), in particular for the minimizer pp. Note that 
/3 = 00 is allowed, c.f. Eq. (|2.22|) . and we can also let p = pp vary for /3 — > 00, 
because II/O/3II5/3 is uniformly bounded by Lemma 17!^ (ii) . Altogether, by 1)2. 13p . 

, P^^(p,T,B,Z) , , 

as Z ^ 00, uniformly in B as long as B/Z^ 0. 
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5.2 Lower bound 

We use the Peierls-Bogoliubov inequality ([HI, (2.1.7), (2.1.8)): If A, B and A + B 
are self-adjoint operators and F{A) = Intre""^/"^, then 

F{A + B)>F{A)-{B)a (5.10) 

where 

We use this inequality with A + B = Hz.b — /uA^ and A = Hz,B,p — IJ-^ — D{p, p). 
Then B is, apart from the constant D{p,p), the second quantization of 

N N 

J^|xi-x,[-i- J];p*[xi|-'- (5.12) 

i<j i=l 

In terms of the creation and annihilation operators al{x.) and as(x) (s= spin com- 
ponent) this can be written 



-B = - J J a*(x)a*/(x')|x — x'l ^as'{x.')as{x) dxdx.' 

s,s' 

-/ / <{^)as{x)\^-y\-^dxdy + D{p,p) (5.13) 

s,s' 

We thus get 

P^^{p,B,Z,T) yTlnUve-^^^'^'"-"^^/^) +D{p,p)-T{B)a (5.14) 



We want to show that we can choose p such that {B)a < 0. 

Since A is the second quantization of a one-particle operator, all expectation 
values of products of creation and annihilation operators in the state {■)a can be 
written in terms of two point correlations by using Wick's theorem. The expectation 
value {B)a involves terms of the form 

(a*(x)a*,(x')a^'(x')as(x))A, (5.15) 

which by can be written 

(a*(x)as(x))A(a*,(x')as'(x'))A - (a*(x)as'(x'))A(a*,(x')as(x))A 

= (a:(x)a,(x))A(a:,(x')a,,(x'))A - |(a:,(x')a,(x))Ap. (5.16) 

If we define 

p(x) = j;(a:(x)a,(x))A (5.17) 
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then, using (|5.16p . we can write 

s,s' 

< D{p-p,p-p). (5.18) 



We now choose p = p, i.e., we choose p = p™^ where p™^ is the solution of the Hartree 
equation H3.4() . This gives 

P^{p,Z,B,V)>P'^^{p;Z,B,T). (5.19) 

By Corollary Ol P'^^ip; Z, B,T)/P^'^^ {p; Z, B,T) converges to 1 in the limit con- 
sidered. Thus 

,^JWVB^ (5.20) 

P^^P{p,T,B,Z) - ^ ' 

and the proof of Theorem 1.1 is complete. 



6 Conclusions 

Using magnetic coherent states we have proved a semiclassical limit theorem for a 
mean field quantum mechanical pressure functional and applied it to derive MFT 
theory at nonzero temperatures as a limit of quantum statistical mechanics in a 
certain parameter range. Our result concerns the grand canonical partition function 
while the corresponding result for the free energy is left as an open problem. Other 
interesting questions not tackled here concern the extension of the asymptotic ground 
state classification to temperature and field strength regions beyond those of the 
present analysis, as well as the thermodynamic limit, where the number of nuclei 
tends to infinity. 
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